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Abstract. We have carried out multi-epoch, time-series differential /-band photometry of a large sample of objects 
in the south-east region of the young (~3 Myr) , nearby (~350 pc) a Orionis open cluster. A field of ~1000 arcmin^ 
was monitored during four nights over a period of two years. Using this dataset, we have studied the photomet- 
ric variability of twenty-eight brown dwarf cluster candidates with masses ranging from the stellar-substellar 
boundary down to the planetary- mass domain. We have found that about 50% of the sample show photometric 
variability on timescales from less than one hour to several days and years. The amplitudes of the /-band light 
curves range from less than 0.01 up to ~0.4 magnitudes. A correlation between the near-infrared excess in the Ks 
band, strong Hq emission and large-amplitude photometric variation is observed. We briefly discuss how these 
results may fit the different scenarios proposed to explain the variability of cool and ultracool dwarfs (i.e. mag- 
netic spots, patchy obscuration by dust clouds, surrounding accretion discs and binarity). Additionally, we have 
determined tentative rotational periods in the range 3 to 40 h for three objects with masses around GOMjup, and 
the rotational velocity of 14 ± 4kms^^ for one of them. 
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1. Introduction 

The photometric variability of brown dwarfs and very low 
mass stars close to the substellar limit has received sig- 
nificant attention during the past few years. The stud- 
ies, mostly performed in the red part of the optical spec- 
trum, using broad-band (Martin & Zapatero Osorio 1997 
Terndrup et al. 1999; Bailer- Jones & Mundt 1999, 2001 
Martin, Zapatero Osorio & Lehto 2001; Gelino et al. 2002 
Joergens et al. 2003) and narrow-band filters (Tinney 
& ToUey 1999; Clarke, Tinney & Covey 2002; Clarke, 
Oppenheimer & Tinney 2002), indicate that approxi- 
mately one third of the late M- and L-type dwarfs show 
variability with amplitudes from ^10 to '--^SO milimagni- 
tudes (mmag). Several of these objects have been reported 
as periodic variables, with periods in the range from half 
an hour to ten days. Photometry in the near-infrared re- 
gion, where cool dwarfs emit most of their output energy, 
has also been performed in young clusters (Carpenter et al. 
2002; Zapatero Osorio et al. 2003— hereafter ZO03) and 
in the field (Bailer- Jones & Lamm 2003; Enoch, Brown 
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& Burgasser 2003). The detected amplitudes of varia- 
tion are in the range 0.05-0.2 mag, and tentative periods 
of ^1.5 to 3h have been reported. On the other hand, 
time-resolved spectroscopy has emphasized the study of 
dust/temperature- or activity-sensitive spectral features 
such as water vapour, TiO, FeH and CrH bands or Ha 
emission, respectively (Nakajima et al. 2000; Hall 2002; 
Bailer- Jones 2002; Burgasser et al. 2002; Liebert et al. 
2003; Clarke, Tinney & Hodgkin 2003). Among the scenar- 
ios proposed to explain the observed levels of photomet- 
ric and spectroscopic variability, we may list cool and/or 
hot corotating magnetically-induced spots in the atmo- 
spheres, heterogeneous or patchy coverage of photospheric 
clouds of grains of solid condensates, (sub)stellar discs sur- 
rounding the objects, and the presence of very low-mass 
companions in close orbits, which may produce eclipses or 
mass-transfer episodes. 

Photometric variability studies of objects with con- 
firmed membership in star clusters with known astrophys- 
ical properties have advantages to those in the field. Age, 
metallicity and distance are the same for the sample of 
cluster members. Because they lie within a well-defined 
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photometric sequence, different luminosities correspond 
to different masses. Hence, the monitoring of a sample of 
cluster objects usually leads to the exploration of the prop- 
erties of members in a wide range of masses. Additionally, 
wide field cameras can simultaneously observe a rela- 
tively large number of sources under the same sky and 
instrumental conditions. At very young ages (e.g., a few 
Myr), fast-rotating completely convective brown dwarfs 
may generate a magnetic field strong enough to produce 
solar-like surface spots or chromospheric activity. Also, 
accretion discs may not yet have been completely dissi- 
pated. Investigation of the rotation of very young brown 
dwarfs is especially important for understanding crucial 
phenomena related to substellar formation, such as angu- 
lar momentum evolution, the relation between disc and 
central object, and disc evaporation. 

Near-infrared and optical wide-field photometric mon- 
itoring of stars in the Orion region has already been 
performed, especially in the area surrounding the Orion 
Nebula Cluster (Herbst et al. 2000, 2002; Carpenter, 
Hillenbrand & Skrutskie 2001; RebuU 2001 and references 
therein). In particular, in the a Orionis region, evidence 
of variability has been found in a large variety of ob- 
jects, from the massive magnetic helium-strong B2Vp star 
in the central OB multiple stellar system, through solar- 
like stars, to objects below the stellar-substellar threshold 
(e.g. Fedorovich 1960; Landstreet & Borra 1978; Bailer- 
Jones & Mundt 2001). 

In this paper we present differential photometry of 
32 low-mass member candidates of the young a Orionis 
open cluster, which belongs to the Ori OBlb association. 
Walraven photometry and hipparcos measurements of 
the distance to the brightest stars of the association point 
to a distance of ~350pc (Brown, de Geus & de Zeeuw 
1994; Ferryman et al. 1997; de Zeeuw et al. 1999). Optical 
extinction towards the a Orionis stellar system and the 
surrounding region seems to be small (Ay <lmag — Lee 
1968; Bejar et al. 2001; Oliveira et al. 2002). The age of 
the cluster is estimated between 1 and 8 Myr, with 3 Myr 
being the most likely value (Zapatero Osorio et al. 2002a; 
Oliveira et al. 2002). 

2. Observations and analysis 

The observations were performed at Observatorio del 
Roque de Los Muchachos, on La Falma (Canary Islands). 
We used the four-chip Wide Field Camera (WFC) mosaic, 
mounted at the prime focus of the 2.5 m Isaac Newton 
Telescope. The detector consists of four thinned EEV4280 
2k X 4k CCDs with a pixel size of 13.5 /im, corresponding 
to a pixel scale of 0.333". 

Monitoring of a ^^0. 29 degree^ area south-east of the 
centre of the a Orionis cluster was carried out during 
two epochs: three consecutive nights, from 2000 December 
30 through 2001 January 1 (hereafter, WFCOO), and on 
2003 January 8 (WFC03). Several Landolt standard stars 
(Landolt 1992) were observed during 2003 January 8, 
which allowed us to calibrate the WFC03 / magnitude. 



Table 1. Log of the observations 



Date 


No. 
exposures 


Airmass 


FWHM 

n 


Exp. time 


30 Dec 2000 

31 Dec 2000 
01 Jan 2001 


5 

10 
6 


1.23-1.70 
1.17-1.50 
1.17-1.58 


1.4-1.6 
0.9-1.4 
1.2-1.8 


1500 
1500 
1500 


08 Jan 2003 


12 


1.17-1.81 


1.0-2.0 


1200 



Thirty-three different exposures, each with completeness 
magnitude / ^^21. 5, resulted from these observations. We 
used the RGO / filter, which resembles the Johnson- 
Cousins-Kron /. Although the nights were clear, seeing 
conditions varied slightly during the runs. In the first 
epoch, a longer exposure time was used, and the counts of 
several bright objects lay in the non-linear regime when 
the seeing improved. In TableHis a log of the observations 
with other useful information. We give the date, number of 
exposures taken each observing night, effective mean air- 
mass, average full width half maximum (FWHM) in arcsec 
and exposure time in seconds. Above the horizontal line 
is shown the WFCOO dataset, and below is the WFC03 
dataset. The WFCOO data have already been used, to- 
gether with near-infrared observations, to study the pho- 
tometric variability of the young, low-mass brown dwarf 
S Ori 45 (ZO03). 

The reduction of the images was performed within 
the IRAF^ environment, using standard packages (imred- 
CCDRED). Raw images were bias-subtracted and flat- 
fielded before alignment. A dither pattern, with displace- 
ments of several times the seeing, was used in order to 
subtract the contribution of the fringing appropriately. To 
accomplish this, a flat-field image for each night was con- 
structed using bias-subtracted scientific images. 

Alignments of the WFCOO and WFC03 images were 
done separately, because of the presence of a 0.6-arcmin 
shift of pointing centres between epochs. Residual fringing 
was below the level of 0.8 % over the CCDs. A study of the 
average FWHM and sky background of each image was 
made before calculating the aperture photometry. This 
analysis was done for each detector individually. We used 
the IRAF DAOPHOT package to perform the photometric 
analysis. Aperture photometry was performed using nine 
different apertures per image: 0.50, 0.60, 0.70, 0.80, 0.90, 
1.00, 1.25, 1.50, and 2.00 times the average FWHM. The 
inner radius of the sky annulus was 5.00 times the average 
FWHM, while the width of the annulus was fixed at 8 
pixels. Instrumental magnitudes for 29514 objects were 
obtained. 



IRAF is distributed by National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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2.1. Differential photometry 

The data analysis (selection of sources, calculation of 
the reference magnitudes, statistical analysis, obtaining 
of light curves, preliminary time-series analysis) was done 
mainly with the matlab^ interactive system. 

We required the mean of the iraf magnitude er- 
rors from sources of the whole sample to be lower than 
0.25 mag for each of the nine apertures. In order to avoid 
source confusion, especially when faint objects are close to 
bright ones or in crowded fields, objects with standard de- 
viations of the centroid coordinates larger than 1.5 pixels 
(~0.5") were rejected. Taking into account these require- 
ments, sources with counts close to saturation, very faint 
objects that appear close to very bright stars, cosmic rays 
and chip defects were mostly discarded. 

Instrumental magnitudes were transformed into differ- 
ential magnitudes, /diff, by subtracting a reference value 
(i.e. a zero-point magnitude level). This reference magni- 
tude for each CCD and exposure was calculated using a 
large number of bright, non- variable stars located over the 
entire monitored area. They were chosen using an itera- 
tive procedure: only low-photometric-noise sources with 
standard deviation cr(/diff) < 0.020 mag (in the WFCOO 
dataset) and < 0.015 mag (WFC03), and with magni- 
tudes in the range 17.5 mag < I < 21.0 mag, were selected 
from an initial sample of reference star candidates. Each 
CCD of the WFCOO and WFC03 datasets was treated 
separately. The final sample of reference objects (~100 
per chip) was the intersection between the WFCOO and 
WFC03 reference samples. We note that a photometri- 
cally stable brown-dwarf target, when compared to other 
objects of similar magnitude, may also be considered a 
reference "star" if it matches the criteria given above. 

Reference magnitudes were calculated following Bailer- 
Jones & Mundt (2001). The method is based in a mean of 
fluxes of the reference stars, where the brightest sources, 
with the lowest Poissonian photon-count errors, dominate 
the final reference magnitude. However, a simple mean of 
the magnitudes of the reference stars was also worked out, 
to check for any significant difference. Because of the large 
size of the sample of reference stars, both zero-point mag- 
nitudes were similar. As the standard deviation of the dif- 
ference between the mean of fluxes and the mean of mag- 
nitudes was lower than 1 mmag, we used only the mean- 
of- fluxes reference level in the subsequent calculations. 

We studied the dependence of the reference magni- 
tude for each CCD on time, airmass and seeing. No 
trends were detected, so we did not applied either airmass- 
or seeing-dependent high-order corrections. Computations 
were made for each Wide Field Camera chip separately, 
as each CCD has different characteristics (i.e. gain, elec- 
tronic noise and offsets between instrumental magnitudes 
and apparent magnitudes). Also, for each chip the average 
FWHM is different at a given time because of the different 
distances to the centre of the focal plane. 





^ MATLAB is a high-performance language for technical com- 
puting developed by The MathWorks, Inc. 
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Right Ascension (h) 

Fig. 1. Spatial location of the 32 brown-dwarf member 
candidates in the a Orionis cluster studied in this article, 
shown with filled squares. Stars brighter than / — 15.0 mag 
that presumably belong to the cluster are represented with 
open circles. The circle radii are inversely proportional to 
the magnitude. Rectangles denote the four CCD chips of 
the Wide Field Camera. 



We performed the previous analysis in the WFCOO 
dataset for each of the nine apertures (0.5 to 2.0 x 
FWHM) and chose the 0.70 x FWHM aperture, which 
both minimized the photon contribution of nearby con- 
taminant sources and maximized the signal-to-noise ratio 
and, as a result, provided the largest number of reference 
stars. 

An additional offset was added to the Idis values of 
each object in order to transform them into apparent mag- 
nitudes, calibrated using the Landolt standard stars ob- 
served during the WFC03 run. As a result, for each object 
the mean WFC03 /diff, denoted by /qs, is equal to the real 
/-band magnitude. 

Finally, multi-epoch calibrated /-band light curves 
were obtained for more than 8000 objects that fulfilled 
all the criteria given above. 

3. Substellar candidates 

Of the 32 (7 Orionis brown dwarf candidates studied here, 
19 were photometrically discovered at optical wavelengths 
by Bejar et al. (1999, 2004). Their near-infrared colours 
are consistent with membership of the cluster. In addition, 
many of them have been confirmed spectroscopically. The 
remaining 13 objects were selected from a combined IJ 
survey conducted with the Wide Field Camera (/) at the 
Isaac Newton Telescope and the isaac instrument (J) at 
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Table 2. Photometric and spectroscopic data for the 32 cr Orionis cluster members considered in our study 



Name 


I 


I -J 


Sp. type 


pEW(Li i) (A) 




pEW(Hq) (A) 












BZOR 


zo 


ByN M 


S Ori J054000.2-025159 


16.21±0.04 


2.15±0.05 


- 


- 


- 


- 


- - 


S Ori J053902.1-023501 


16.41±0.03 


1.62±0.05 


- 


- 


- 


- 


- 


S Ori 16" 


16.62±0.03 


1.94±0.04 


- 


- 


- 


- 


- 


S Ori J053847.2-025756 


16.75±0.02 


2.23±0.03 


- 


- 


- 


- 


- 


S Ori J053829.0-024847'' 


16.81±0.02 


1.80±0.06 


M6.0±0.5 


- 


- 


- 


< 5 


S Ori J053954.3-023719" 


16.82±0.03 


1.94±0.04 


- 


- 


- 


- 


- 


S Ori J053825.4-02424r 


16.86±0.02 


1.73±0.03 


- 


- 


- 


- 


- 


S Ori 25" 


16.93±0.03 


2.16±0.03 


M7.5±0.5 


0.6 


45.0±1.0 


- 


42±8 44±1 


S Ori J053826.1-02404r 


16.93±0.02 


1.95±0.04 


M8.0±0.5 


- 


- 


- 


4±2 


S Ori 27" 


17.03±0.03 


2.15±0.03 


M7.0±0.5 


0.74±0.09 


6.1±1.0 


5.7±0.5 


5±2 


S Ori J053922.2-024552 


17.05±0.03 


1.65±0.05 


- 


- 


- 


- 


- 


S Ori J054014.0-023127 


17.07±0.03 


2.60±0.04 


- 


- 


- 


- 


- 


S Ori 28" 


17.14±0.03 


1.81±0.03 


- 


- 


- 


- 


- 


S Ori 31" 


17.23±0.02 


1.87±0.07 


M7.0±0.5 


- 


- 


- 


2.5±0.9 


S Ori 30" 


17.29±0.03 


1.98±0.03 


M6.0±0.5 


- 


- 


- 


16±6 


S Ori 32" 


17.34±0.03 


1.85±0.03 


- 


- 


- 


- 


- 


S Ori J054004.5-023642" 


17.51±0.03 


2.02±0.04 


- 


- 


- 


- 


- 


S Ori 36" 


17.71±0.03 


2.20±0.04 


- 


- 


- 


- 


- 


S Ori J053918.1-025257 


18.64±0.04 


2.51±0.04 


- 


- 


- 


- 


- 


S Ori 42" 


19.04±0.03 


2.48±0.03 


M7.5±0.5 


- 


- 


- 


89±12 


S Ori 45" 


19.49±0.02 


2.57±0.02 


M8.5±0.5 


2.4±1.0 


60.0±1.0 


33±9 


26±15 


S Ori J053929.4-024636 


19.73±0.03 


2.55±0.03 


- 


- 


- 


- 


- 


S Ori 71" 


20.02±0.03 


2.69±0.04 


L0.0±0.5 


- 


- 


- 


700±80 


S Ori J053849.5-024934 


20.08±0.02 


2.80±0.04 


- 


- 


- 


- 


- 


S Ori 51" 


20.23±0.03 


2.85±0.07 


M9.0±0.5 


- 


- 


- 


25: 


S Ori 50" 


20.48±0.03 


2.88±0.04 


M9.0±0.5 


- 


- 


- 


< 10 


S Ori 47" 


20.49±0.02 


3.16±0.03 


L1.5±1.0 


4.3±0.5 


- 


< 6 


25: 


S Ori J053944.5-025959 


20.74±0.04 


2.72±0.05 


- 


- 


- 


- 


- 


S Ori 53" 


20.96±0.02 


2.87±0.03 


M9.0±0.5 


- 


- 


- 


< 10 


S Ori J054007.0-023604 


21.16±0.03 


2.84±0.04 


- 


- 


- 


- 


- 


S Ori J053956.8-025315 


21.25±0.04 


2.97±0.05 


- 


- 


- 


- 


- 


S Ori J053858.6-025228 


22.17±0.05 


3.54±0.06 


- 


- 


- 


- 


- 



References: Bejar, Zapatero Osorio & Rebolo 1999 and Bejar et al. 2004 (BZOR); Zapatero Osorio et al. 1999, 2002a (ZO); 
Barrado y Navascues et al. 2001, 2002, 2003 (ByN); MuzeroUe et al. 2003 (M). 



the Very Large Telescope Antu/UT 1. Details of this sur- 
vey will be given in a forthcoming paper. The 32 selected 
targets for the present photometric analysis have / mag- 
nitudes between 16.0 and 22.5 mag. Neither brighter nor 
fainter objects yielded reliable differential photometry ei- 
ther because of non-linearity or excessively faint signals in 
each exposure, respectively. In Fig.Q] we show the spatial 
location of the targets, stars brighter than / — 15.0 mag 
that presumably belong to the a Orionis cluster (based on 
spectral type and/or proper motion determinations) and 
the position on the sky of the WFC mosaic. We provide in 
TableElthe complete sample, together with relevant spec- 
troscopic information from the literature (e.g., spectral 
types and pseudo-equivalent widths of lithium absorption 
at 6708 A and Ha emission). 



The spectral types of the targets range from ~M5.5 
to '^L2. According to theoretical models, at the age of 



the cluster (~3 Myr), these roughly correspond to masses^ 
between ~80 and ~8Mjup. The upper hmit is shghtly 
larger than the theoretical minimum mass for hydrogen fu- 
sion, the borderline between stars and brown dwarfs. The 
lower limit is below the deuterium-burning mass thresh- 
old, suggested as the frontier between brown dwarfs and 
planetary-mass objects (e.g. Saumon et al. 1996). These 
mass estimates may be affected by uncertainties of the 
evolutionary models at very young ages (see the discus- 
sion by Baraffe et al. 2002, 2003). The spectrophotometric 
sequence of the a Orionis cluster in the substellar domain 
is known from previous studies (references in footnote of 
Table 121 and Bejar et al. 2001). A colour-magnitude dia- 
gram I vs. I— J, such as the one shown in Fig. 13 provides a 
reliable way of determining the membership of candidates 
without spectroscopic data. 



Conversion between solar and Jupiter masses: 
1047.56 ± 0.08 Mjup. 



IMp 
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Fig. 2. /vs. I — J colour-magnitude diagram of 
the 32 brown-dwarf member candidates under study. 
Spectroscopically confirmed substellar cluster members 
are shown with filled squares, while open squares indicate 
objects with no available spectroscopic data. The solid 
and dashed lines are 3 Myr colour-magnitude isochrones 
for substellar objects, based on models by the Lyon group 
(BarafFe et al. 2002 and Chabrier et al. 2000, respectively). 
Approximate masses are provided to the right. 



3.1. Ha emission 

According to the standard scenario of classical T Tauri 
stars, broad Ha emission lines are caused by the accretion 
flow of free-falling gas. The infalling gas, in this case, is 
accreted during the early stages of star formation on to 
the disc or at the base of high-latitude magnetic accretion 
columns on the stellar surface (Konigl 1991, and references 
therein). Several authors have suggested the extrapolation 
of this picture to objects below the stellar/substellar mass 
limit. 

The great majority of the spectroscopically confirmed 
brown dwarfs show Ha emission (see Table EJ, which may 
be related to chromospheric activity or to the presence of a 
circum(sub)stellar disc. The asymmetric, extraordinarily 
strong Ha emission of the cool substellar object S Ori 
71 should be noted. Barrado y Navascues et al. (2002) 
estimated its mass to be between 22 and 13Mjup. They 
measured a pseudo-equivalent width of 700 ± 80 A, which 
is among the strongest Ha emissions so far seen in a young 
late-type dwarf. S Ori 25, S Ori 42 and S Ori 45 also 
show strong Ha emission, which is ~10-20 times greater 
than the average emission of other a Orionis sources of 
similar spectral type. According to Barrado y Navascues & 
Martin (.2003) , such strong emission is consistent with disc 



Fig. 3. I — J vs. I — Ks colour-colour diagram. All 
the brown-dwarf candidates studied with 2MASS JHKg 
photometry data with precision better than 0.1 mag are 
shown. The isochrones are described in the caption of Fig. 

m 



accretion phenomena. Other cluster objects show weak to 
moderate Ha emission lines. 

Lamm et al. (2004) have recently detected a strong cor- 
relation between photometric variability in periodic vari- 
ables and Ha emission in a sample of pre-main sequence 
stars in the young (2-4 Myr) open cluster NGC 2264. 
However, no evidences for a correlation for the irregular 
variables have been found. 



3.2. Near-infrared excess 

We have used the 2MASS All-Sky Catalog of Point 
Sources (Cutri et al. 2003) to study the presence of any 
near-infrared excesses in our targets that could be an in- 
dication of warm surrounding discs. In Fig. |21 is shown 
the I — J vs. I — Kg diagram for 17 of the brown- 
dwarf candidates studied, with 2MASS JHKs photom- 
etry more precise than 0.1 mag (i.e. targets brighter than 
/ =17.5 mag, approximately). The solid and dashed lines 
are the 3 Myr colour-colour isochrones for substellar ob- 
jects, based on models by the Lyon group (BarafFe et al. 
2002 and Chabrier et al. 2000, respectively). AU the sub- 
stellar candidates except two lie in the dashed isochrone 
within the error bars. The outliers are S Ori 27 and S 
Ori J053825.4-024241. S Ori 27 is an M6.5 ± 0.5 brown 
dwarf with an estimated mass of ~40Mjup. According 
to its / - J colour, S Ori J053825.4-024241 may be an 
^M6-type brown dwarf for which optial spectroscopy is 
not available. While S Ori 27 seems to move towards bluer 
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Fig. 4. J — H vs. H ~ Ks colour-colour diagram. Small 
filled squares with dotted error bars are for spectroscop- 
ically confirmed objects of the a Orionis cluster. Open 
squares with solid error bars denote the objects studied in 
this paper. Spectroscopically confirmed objects are plot- 
ted as filled squares. The dashed line is the 3 Myr colour- 
colour isochrone for substellar objects, based on Chabrier 
et al. (2000) models. Late-K and early-M stars are located 
in the upper left corner of the figure, while the mid- and 
late-M dwarfs are close to the center. 



colours, S Ori J053825. 4— 024241 shows a clear reddening 
that suggests an excess of infrared emission at 2.2 /im or, 
alternatively, it may be affected by interstellar extinction. 
These objects do not lie close to the cluster sequence in 
the J — _ff vs. H — Ks diagram, as shown in Fig. ^ where, 
again, only 2MASS JHR^ data more precise than 0.1 mag 
have been plotted. 

Recently, Oliveira et al. (2003) have obtained L' pho- 
tometry of a Orionis stars with late K and early M 
spectral types, concluding that about 50% of the stellar 
cluster population show a considerable infrared excess at 
3.8 /im. About half of the a Orionis stars harbour dusty 
envelopes that may occur in the form of discs. We note 
that the frequency of strong Ha emissions in a Orionis 
stars and brown-dwarf members (Zapatero Osorio et al. 
2002a, 2002b; Barrado y Navascues 2003; MuzeroUe et al. 
2003) is not as high as the frequency of L'-band infrared 
excesses. Furthermore, the rate of K-hand infrared ex- 
cesses is rather low as compared to the i'-band (Oliveira 
et al. 2002). Mid-infrared excesses are better disc indica- 
tors. 

Evidence for accretion discs, dusty envelopes and 
the existence of a T Tauri-like phase in young brown 
dwarfs has been claimed for other similarly young clusters 



and star-forming regions, such as p Ophiuchi (Wilking, 
Greene & Meyer 1999; Testi et al. 2002; Natta et al. 
2002), the Trapezium (Muench et al. 2001), Chamaeleon 
I (Natta & Testi 2001; Persi et al. 2001; Jayawardhana, 
Mohanty & Basri 2002), R Coronae Australis (Fernandez 
& Comeron 2001, Barrado y Navascues, Mohanty & 
Jayawardhana 2003), IC 348 (Jayawardhana, Mohanty & 
Basri 2003), Lupus 3 (Comeron et al. 2003), and TW 
Hydrae (Mohanty, Jayawardhana & Barrado y Navascues 
2003). 



4. Criteria of variability 

The scatter of the data points in a light curve repre- 
sents a measure of the photometric variability, or at least 
sets an upper limit to the variability. The standard de- 
viation, (t(/), of the light curve is the most used vari- 
ability indicator, although other parameters have been 
also proposed: the peak-to-peak amplitude, the average 
of the absolute relative magnitudes, the "range" (the dif- 
ference in median values of the upper 15 % and the lower 
15 % of measurements), variability and median statistics, 
as the x^ or the fj, and other indices or tests as the 
Kolmogorov-Smirnov test. Some of the parameters need 
posterior Monte Carlo simulations to quantify the like- 
lihood of variability. Discussions about some variability 
indicators can be found in Stetson (1996), Bailer-Jones 
& Mundt (1999), Herbst et al. (2002) and Enoch et al. 
(2003). Here, we will define a few criteria of variability ac- 
cording to different timescales. We will compare the pho- 
tometry of our a Orionis sample with the data of field 
sources of similar brightness. For consistency purposes, 
all of the light curves are treated in the same manner. 

We define a{I) as the standard deviation of the light 
curves for each night. The a{I) vs. / diagram shown in 
Fig. 13 illustrates the strong dependency of the photomet- 
ric standard deviation on the target apparent magnitude. 
The lower envelope of a{I) is roughly fiat at the level of 
5 mmag in the / ^17.0-19.5mag range, while it increases 
exponentially above '^20.5mag because of the Poissonian 
photon-count error. On the other hand, the scatter also 
increases below ^16.5 mag because of saturation and non- 
linear effects in bright sources. We note that to compute 
the (t(/), three / values in the WFCOO epoch and one 
value in the WFC03 were rejected for all the objects, as 
the measured point spread functions (PSF) in the corre- 
sponding images were far from the average PSF for each 
epoch. 

The three brightest brown dwarfs in the sample of 32 
a Orionis objects (/ <16.7mag) will not be considered in 
the discussion of Sect. 5 because a significant fraction of 
the light curve points was recorded with counts in the non- 
linear regime of the detector. Likewise, the light curve of 
S Ori J053858.6-025228, a substellar object fainter than 
/ = 22 mag, was discarded because of the large Poissonian 
error. For the remaining 28 brown dwarfs, what we call the 
final sample, the photometry is reliable. 
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Fig. 5. Light curves of the fourteen brightest objects in the final sample. 



The light curves of our final sample are displayed in 
Figures and El For comparison purposes, we have used a 
common display: vertical height has been fixed to 0.5 mag 
and the temporal scale is constant along the horizontal 
axis. The Modified Juhan Date (MJD = JD -2400000.5) 
is indicated below, and the two epochs are shown. We have 
not taken into account the heliocentric correction, as the 
maximum difference between JD and HJD (Heliocentric 
Julian Date) is of the order of only a few seconds over a 
period of three nights. The reference star light curves are 
similar to that of S Ori 36. 



4.1. Short timescales 

In order to find photometric variability on very short 
timescales (in a series of consecutive exposures, i.e., min- 
utes to hours), we compared the a{I) of the hght curves 
of the targets with those of field stars of similar magni- 
tude. For each target we defined a short-term variability 
parameter, zust, given by the following equation: 

N- ... 

"^""^^ N+ + N-' ^^> 

where the subscript "ST" stands for short-term, and N'^ 
and N~ are the number of sources in the same magni- 
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Fig. 6. Light curves of the fourteen faintest objects in the final sample. 



tude bin with cr(/) larger and lower, respectively, than 
the a{I) of the object under consideration. High values of 
njsT imply more likely variability. The width of the bin 
was fixed to 0.2 magnitudes, centred on the / magnitude 
of the object. Typically about one hundred objects with 
similar brightness and, hence, with similar photometric er- 
ror, were considered in the comparison. Such a narrow bin 
was selected to minimize the error contribution from faint 
sources and sources in the non-linear domain, which show 
a larger photometric error that could mask the intrinsic 
variability of the target. 



We will impose a cut at ccjst — 0.95 to separate vari- 
able from non-variable sources. Targets with tust ^ 0.95 
are considered to be very likely short-term variable brown 
dwarfs in our sample. This criterion demands that for a 
target brown dwarf to be detected as variable it must 
be significantly more variable than field objects in the 
same magnitude bin. Statistically, the light curves of the 
variable sources display photometric amplitudes that are 
among the 5 % largest amplitudes for their magnitude 
bins. This is as if there is only a 5% chance that a brown 
dwarf which is not truly variable is incorrectly classified as 
variable (i.e. 5% chance of a false positive detection). We 
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Fig. 7. Standard deviation of the differential light curves 
versus the / magnitudes for each target (small dots) in the 
second WFCOO night. Only objects with a{I) < 0.15 mag 
are shown. Mean photometric iraf errors are below 5,15 
and lOOmmag for objects brighter than 7=18, 20 and 22, 
respectively. Filled squares denote the brown-dwarf can- 
didates. Short-term variable brown dwarfs (on the basis 
of the definition of wst) are labelled. 

note that the distribution of tngT is rather flat between 
and 1 within a given magnitude bin. However, the distri- 
bution of C7{I) is asymmetric: most of the data points fall 
close to the mean photometric iraf errors, while only a 
small percentage lies far from this location (see Fig.d). 

4.2. Intermediate timescales 

To detect variability in a target from one night to another 
(hereafter mid-term variability, i.e., scale of a few days), 
we compared the averaged magnitudes (of both targets 
and field sources), li, computed for each night, i. We hst 
the mean magnitudes of our sample in Table |3 The com- 
parison is performed by obtaining the absolute value of 
li — Ij per object, where i and j represent two different 
nights. In order to avoid the bias introduced by short-term 
photometric variability, we have divided the differences by 
the standard error of the mean magnitudes. Hence, we 
have defined a dimensionless absolute magnitude differ- 
ence, 5{I)ij, as: 



6{I\ 



\h - h 



<r{uy 



Ahr 



(2) 



where n stands for the number of photometric data points 
per night. By using mean /-band magnitudes per night we 
are increasing the sensitivity to mid-term variability. 



Fig. 8. Dimensionless absolute differences of mean mag- 
nitudes between the second and third WFCOO nights vs. 
apparent / magnitudes. Field objects are plotted as small 
dots, while filled squares denote the brown-dwarf can- 
didates. Names of several mid-term variable targets are 
given. 



We have selected likely mid-term variables in our sam- 
ple in a similar manner as we picked short-term variable 
brown dwarfs. We have used the definition of equation ^ 
which we now name tumt, where the subscript stands for 
"mid-term" . 7V+ and N~ represent, in this case, the num- 
ber of sources in the same magnitude bin with 5{I)ij larger 
and lower, respectively, than the S{I)ij of the object un- 
der consideration. Figure |S1 shows the dimensionless dif- 
ferences against apparent magnitudes for the second and 
third nights of the WFCOO epoch. Similar diagrams are 
obtained when considering other possible combinations, 
e.g. first and second nights, second and third nights. We 
note that mid-term variability can be investigated only 
in the WFCOO data. As in the previous analysis of short- 
term variability, the bin magnitude is chosen to be 0.2 mag 
centred on the / magnitude of each target. According to 
our criterion, likely mid-term variable brown dwarfs have 
wuT > 0.95. 



4.3. Long timescales 

For long-term variability, which takes into account differ- 
ences between the first epoch of observations (WFCOO) 
and the second (WFC03), separated by about 2 years, we 
applied a criterion similar to that of mid-term variabil- 
ity. Now, the only WFC03 night is compared to each of 
the three WFCOO nights using the dimensionless magni- 
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Fig. 9. Dimensionless absolute differences of mean magni- 
tudes between the second WFCOO night and the WFC03 
night vs. apparent / magnitudes. Field objects are plotted 
as small dots, while filled squares denote the brown-dwarf 
candidates. Names of several variable targets are given. 



tude differences, S{I)ij. The parameter that accounts for 
long-term variability is named tjjut- 

Long-term variable a Orionis brown dwarfs in our sam- 
ple have wlt > 0.95. The wlt calculated for the com- 
parison between the first WFCOO night and the WFC03 
epoch is depicted as a function of apparent / magnitude 
in Fig. El Similar diagrams are obtained for the remaining 
two WFCOO nights. We note that according to the adopted 
criterion of variability, all mid-term variables could also be 
identified as long-term variable sources. This is so because 
we have performed comparisons of each WFCOO night to 
the WFC03 epoch. Mid-term variable sources might be 
true long-term variables but we cannot discriminate eas- 
ily with our data. Hence, we have excluded mid-term vari- 
ables from our analysis of long-term variability. 

5. Discussion 

We provide in Table01the variability parameters obtained 
for the three different timescales and for our target sample. 

5.1. Long-term variable brown dwarfs 

Two brown dwarfs, namely S Ori J053847.2-025756 and 
S Ori 42, vary significantly {wlt > 0.95) from the WFCOO 
to the WFC03 epoch. The amplitude of their light curves 
is given in Tabled These objects represent about 7% of 
the brown-dwarf final sample. Because our observations 
are not equally sensitive to /-band long-term variability 



as they are to mid- and short-term photometric variabil- 
ity, this percentage represents a lower limit to the true 
frequency of long-term variable young brown dwarfs. 

The exact temporal scale of the observed /-band long- 
term variability, which may be between several days and 
years, is unknown to us. This kind of variability could 
be associated with the long-term evolution of surface fea- 
tures, such as dust clouds or solar-like spots, or with the 
presence of an accretion disc or a close interactive compan- 
ion. In the case of S Ori J053847.2-025756, spectral data 
are lacking, and its Ha emission is unknown. Given its 
relatively high effective temperature ('^2800K, obtained 
from its I— J colour), the formation of magnetically-driven 
spots or the presence of an accretion disc may be favoured 
over dust condensation (Mohanty et al. 2002; Joergens 
et al. 2003). The brown dwarf S Ori 42 is ~2mag fainter 
and cooler, with an estimated effective temperature at 
around 2500 K. Dust condensation is expected at these 
cool temperatures; however, the notably high Ha emis- 
sion of S Ori 42 (Table[2J) supports the scenarios involving 
mass infalling. 

Among our final list of brown dwarfs, S Ori 71 is the 
strongest Ha emitter (Table EJ. This source is a non- 
variable object according to our criteria. Nevertheless, the 
three ii^lt are quite close to the 0.95 limit, suggesting that 
it may be a real long-term variable with an amplitude of 
~0. 13 mag. S Ori 71 is quite faint, and relatively large pho- 
tometric errors have likely prevented us from detecting a 
clear magnitude variation. The persistent and strong Ha 
emission of S Ori 71 suggests that, for this particular ob- 
ject, the emission is probably linked to the presence of an 
accretion disc or to a close interactive companion (Barrado 
y Navascues et al. 2002). 

5.2. Mid-term variable brown dwarfs 

Six young brown dwarf candidates show significant night- 
to-night /-band variations during the WFCOO run. They 
display at least one of the three tUMT values larger than 
0.95. The amplitude of their light curves is provided in 
Table |S| These six mid-term variable sources represent 
about 21 % of our final sample. Possible causes that lead 
to mid-term photometric variability are: variable obscu- 
ration by circumsubstellar dust, heterogeneous distribu- 
tion of hot or cool spots that cover a large percentage 
of the substellar surface, dust condensates in the atmo- 
sphere, variation in accretion rate and transits of faint 
companions in close orbits. We highlight two cases: S Ori 
J053825.4-024241 and S Ori 25. 

S Ori J053825.4-024241 is the most variable brown 
dwarf candidate in our final sample. The amplitude of the 
light curve, as measured from peak-to-peak, is 0.36 mag. 
As indicated in SectionESl the 2MASS /Cg-band data of S 
Ori J053825. 4— 024241 show significant near-infrared ex- 
cess, which may be indicative of the presence of a sur- 
rounding disc, from which the central object may be ac- 
creting. Hence, the observed /-band variability could be 
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Table 3. Mean magnitude and standard deviation for the light curves of each night" 
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17.15±0.03 


S Or 


31 




17.207±0.016 


17.230±0.017 


17.252±0.007 


17.226±0.019 


S Or 


30 




17.273±0.016 


17.297±0.016 


17.284±0.010 


17.290±0.010 


S Or 


32 




17.33±0.02 


17.300±0.017 


17.356±0.013 


17.34±0.02 


S Or 


J054004.5-^ 


-023642 


17.579±0.012 


17.535±0.013 


17.506±0.007 


17.516±0.020 


S Or 


36 




17.707±0.009 


17.710±0.009 


17.710±0.004 


17.713±0.008 


S Or 


J053918.1- 


-025257 


18.62±0.05 


18.68±0.05 


18.68±0.05 


18.64±0.04 


S Or 


42 




19.064±0.009 


19.142±0.019 


19.088±0.009 


19.04±0.03 


S Or 


45 




19.484±0.013 


16.49±0.03 


19.46±0.03 


19.48±0.03 


S Or 


J053929.4- 


-024636 


19.74±0.04 


19.76±0.03 


19.78±0.04 


19.73±0.02 


S Or 


71 




19.91±0.04 


19.96±0.02 


19.901±0.012 


20.02±0.02 


S Or 


J053849.5- 


-024934 


20.150±0.017 


20.18±0.04 


20.118±0.017 


20.081±0.019 


S Or 


51 




20.22±0.07 


20.24±0.04 


20.23±0.05 


20.23±0.02 


S Or 


50 




20.54±0.09 


20.53±0.03 


20.60±0.04 


20.48±0.02 


S Or 


47 




20.50±0.03 


20.52±0.04 


20.52±0.06 


20.49±0.03 


S Or 


JQ53944.5- 


-025959 


20.78±0.06 


20.75±0.03 


20.75±0.05 


20.74±0.03 


S Or 


53 




20.95±0.06 


21.01±0.03 


20.91±0.06 


20.96±0.04 


S Or 


J054007.0- 


-023604 


21.16±0.06 


21.17±0.07 


21.18±0.08 


21.16±0.07 


S Or 


J053956.8- 


-025315 


21.27±0.04 


21.29±0.05 


21.24±0.03 


21.25±0.07 


S Or 


J053858.6- 


-025228 


21.91±0.20: 


22.23±0.18: 


22.3±0.2: 


22.17±0.11: 



00#1: first observing night of WFCOO; 00#2: second observing niglit of WFCOO; 00#3: third observing night of WFCOO; 
03: observing night of WFC03 



related to episodes of mass accretion or to eclipses caused 
by inner portions of the disc, as those presumably ob- 
served in the pre-main-sequence object KH 15D (Hamilton 
et al. 2001). Herbst et al. (2002) stated that a correla- 
tion between near-infrared excesses and large photomet- 
ric variability in low-mass T Tauri stars exists. S Ori 
J053825. 4— 024241 possibly represents one of the first ex- 
amples of a similar correlation among young brown dwarfs 
(see also those presented by Carpenter et al. (2001) in the 
Chamaeleon I molecular cloud). However, we note that 
the connection between near-infrared excesses and the ex- 
istence of disc accretion is not so straight forward. 

Optical spectra are also available for S Ori 25. They 
show strong Ha emission, with a width of more than 
100 km s^^ at 10% of peak intensity (see references given 
in Table |3l. The Ha emission of S Ori 25 is remarkably 
stable on a timescale of years. In addition, this object 
appears to show a modulation in the light curve with a 
period of ~40h, as we will discuss below. 



The other four mid-term variable brown dwarf can- 
didates are S Ori J054014.0-023127, S Ori 28, S Ori 
32 and S Ori J054004.5-023642. Spectroscopy is not 
available for any of them. However, we note that S Ori 
J054014. 0—023127, for which we determine a photometric 
/-band amplitude of ~0.12mag, appears quite overlumi- 
nous in the colour-magnitude diagram of Fig. |21 

5.3. Short-term variable brown dwarfs 

Three a Orionis brown dwarf candidates are classifed 
as very likely short-term variables according to our 
K7ST criterion. They are S Ori J053825.4-024241, S 
Ori J053826. 1-024041 and S Ori J053918. 1-025257. 
Additionally, although the wst values of S Ori 27 and 
S Ori 28 are below the tost = 0.95 boundary, we have 
measured reliable peaks in the periodograms of their light 
curves that also suggest classifying them as periodic short- 
term variable brown dwarfs (see Sect. 15. 41) . The amplitudes 
of the light curves are given in Table \S\ 
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Table 4. Long-, mid- and short- term variability parameters of the final sample. See Section 4 for the definition of 
each zu 



Name 


tULT 


tJ7LT 


CClt 


zumt 


zumt 


ujmt 


ZUST 


ZUST 


tUST 


tUST 




00#l-03 


00#2-03 


00#3-03 


00#l-00#2 


0()#l-0()#3 


()0#2-()0#3 


00#1 


00#2 


00#3 


03 


S Ori J053847.2-025756 


0.93 


0.97 


0.95 


0.2 


0.5 


0.1 


0.5 


0.2 


0.4 


0.4 


S Ori J053829.0-024847 


0.6 


0.9 


0.8 


0.5 


0.6 


0.1 


0.8 


0.5 


0.9 


0.5 


S Ori J053954.3-023719 


0.6 


0.5 


0.4 


0.0 


0.1 


0.0 


0.5 


0.1 


0.3 


0.8 


S Ori J053825.4-024241 


0.96 


1.00 


0.8 


0.98 


0.99 


0.99 


0.6 


0.98 


0.4 


0.8 


S Ori 25 


0.97 


0.6 


0.9 


0.97 


0.94 


0.94 


0.6 


0.4 


0.6 


0.5 


S Ori J053826.1-024041 


0.4 


0.8 


0.8 


0.8 


0.91 


0.6 


0.94 


0.95 


0.9 


0.6 


S Ori 27 


0.1 


0.8 


0.4 


0.7 


0.4 


0.90 


0.7 


0.7 


0.7 


0.7 


S Ori J053922.2-024552 


0.3 


0.5 


0.6 


0.2 


0.6 


0.6 


0.3 


0.5 


0.2 


0.4 


S Ori J054014.0-023127 


0.2 


0.98 


0.97 


0.96 


0.98 


0.9 


0.4 


0.5 


0.9 


0.8 


S Ori 28 


0.93 


0.6 


0.92 


0.8 


0.3 


0.95 


0.6 


0.5 


0.6 


0.7 


S Ori 31 


0.7 


0.1 


0.6 


0.6 


0.8 


0.8 


0.6 


0.7 


0.2 


0.5 


S Ori 30 


0.6 


0.2 


0.2 


0.6 


0.3 


0.6 


0.5 


0.6 


0.3 


0.2 


S Ori 32 


0.5 


0.8 


0.4 


0.7 


0.7 


0.97 


0.6 


0.7 


0.4 


0.6 


S Ori J054004.5-023642 


0.94 


0.7 


0.3 


0.8 


0.95 


0.90 


0.3 


0.5 


0.2 


0.5 


S Ori 36 


0.2 


0.1 


0.1 


0.1 


0.2 


0.0 


0.2 


0.3 


0.0 


0.1 


S Ori J053918.1-025257 


0.7 


0.7 


0.7 


0.9 


0.94 


0.4 


0.95 


0.96 


0.93 


0.9 


S Ori 42 


0.6 


0.98 


0.8 


0.9 


0.5 


0.8 


0.1 


0.97* 


0.1 


0.5 


S Ori 45 


0.1 


0.0 


0.6 


0.1 


0.6 


0.7 


0.1 


0.7 


0.7 


0.6 


S Ori J053929.4-024636 


0.3 


0.6 


0.8 


0.3 


0.7 


0.6 


0.7 


0.7 


0.7 


0.3 


S Ori 71 


0.9 


0.9 


0.9 


0.6 


0.2 


0.7 


0.5 


0.4 


0.3 


0.2 


S Ori J053849.5-024934 


0.8 


0.94 


0.5 


0.3 


0.5 


0.7 


0.1 


0.6 


0.1 


0.2 


S Ori 51 


0.1 


0.2 


0.0 


0.2 


0.1 


0.1 


0.8 


0.8 


0.5 


0.2 


S Ori 50 


0.7 


0.6 


0.8 


0.3 


0.7 


0.8 


0.90 


0.4 


0.3 


0.1 


S Ori 47 


0.2 


0.6 


0.4 


0.3 


0.3 


0.0 


0.1 


0.7 


0.6 


0.1 


S Ori J053944.5-025959 


0.4 


0.1 


0.1 


0.4 


0.4 


0.1 


0.6 


0.4 


0.4 


0.1 


S Ori 53 


0.2 


0.8 


0.4 


0.6 


0.4 


0.7 


0.4 


0.1 


0.4 


0.3 


S Ori J054007.0-023604 


0.0 


0.3 


0.1 


0.2 


0.2 


0.0 


0.3 


0.7 


0.5 


0.7 


S Ori J053956.8-025315 


0.2 


0.4 


0.1 


0.3 


0.3 


0.3 


0.1 


0.2 


0.2 


0.5 



The S Ori 42 light curve was affected by a crossing artificial satellite path. 



The strong intra-night photometric variations 
observed in S Ori J053826. 1-024041 and S Ori 
J053918. 1-025257 lasts for three or more days. Both 
sources are apparently stochastic photometric variables 
with peak-to-peak amplitudes of ~0.10-0.15mag on 
timescales as short as the time interval between consec- 
utive exposures (~20min). We cannot classify these two 
objects as mid-term variables according to our criteria. 
However, S Ori J053825. 4-024241, which has an infrared 
excess emission at 2.2 /xm, complies with our criteria of 
mid-term variability. 

In principle, short-term variability is more likely to be 
related to fast variations on the surface of an object than 
in a disc, where the relevant timescales are likely to be of 
the order of the revolution period of its particles. However, 
the rotational modulation scenario requires very fast spin 
periods, as short as one or two hours. In contrast, intense 
vertical motions of convcctive cells (bubbling) caused by 
turbulent convection could modify the photospheric cov- 
erage of dust clouds or dark magnetic spots on very short 
timescales. This is supported by the greater importance 
of the internal heat source during the early stages of evo- 



lution of substellar objects, which pushes warm material 
from the interior to the photosphere. Mass infalling from a 
disc and formation of a hot spot could also explain the sud- 
den high-amplitude brightening (by more than 200 mniag) 
of the unsteady variable S Ori J053825.4-024241 during 
the second WFCOO night. 



5.4. A search for periodicities 

Time-series analysis has been performed over our final 
sample of brown dwarfs, i.e. 56 light curves (28 brown- 
dwarf candidates x 2 epochs). We have used period, a 
time-series analysis package (developed by Vik Dhillon 
on La Palma) within the STARLINK software collection, 
which includes the Lomb-Scargle (Lomb 1976; Scargle 
1982) and the CLEAN (Roberts, Lehar & Dreher 1987) 
algorithms. Frequency search limits were the Nyquist crit- 
ical frequency and the inverse of the maximum temporal 
coverage. Five iterations and a gain of 0.1 were used dur- 
ing the CLEAN runs. The periodograms of three targets 
(namely, S Ori 25, S Ori 27 and S Ori 28) show significant 
peaks, which are at least five times higher than the typical 
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Fig. 10. WFCOO (solid line) and WFC03 (dashed line) 
light curve periodograms of S Ori 25, S Ori 27 and S Ori 28 
(from top to bottom) . Insets show the phased light curves 
associated with each peak of the periodogram (WFCOO, 
WFC03 and WFC03, respectively). Note that the scale of 
the CLEANed power and / magnitude and the epoch of S 
Ori 25 are different from S Ori 27 and 28. Further details 
are given in the text. 

noise of the periodograms of non-variable field objects of 
similar brightness. 

S Ori 25 is a mid-term variable brown dwarf with a 
light curve that seems to rise and fall in the three-night 
WFCOO run. The CLEANed periodogram in the top box of 
the FigurelTUIshows a very powerful peak at 40 ± 8h (0.61 
± 0.13 day^^) and its harmonics. The peak also appears 
in the Lomb-Scargle periodogram. The WFC03 data are 
not useful for checking this detection, as the coverage was 
only over ^5h. We will adopt a radius of 0.45 ± 0.05 i?©** 
for S Ori 25. The uncertainty in the radius determination 
is estimated from the age interval 1-8 Myr and the evo- 
lutionary tracks of Chabrier et al. (2002). The projected 
rotational velocity, usini, of 9.4 ± l.Okms"^ (measured 
by MuzeroUe et al. 2003) gives a maximum rotational pe- 
riod of P/sini = 58 ± 9h for S Ori 25, which is com- 
patible with our results. Furthermore, using the measured 
rotational period, we can constrain both the inclination 
angle {i — 'iG'^]^ degrees) and the rotational velocity (14 
± 4kms^^). The 40 h period is greater by far than those 
in the literature for field brown dwarfs (except for the 
238 h period detected in 2MASSW J1300425-I-191235 by 
Gelino et al. 2002), but similar to those recently discov- 
ered by Joergens et al. (2003) in the young Chamaeleon 
I star-forming region. We note that while the duration of 



IRp 



6.9599 X 10** m. 



the periodic modulation of S Ori 25 appears to be gen- 
erally longer than that of field brown dwarfs, its wsini 
measurement is similar (Mohanty & Basri 2003; Bailer- 
Jones 2004). In addition, the strong Ha emission in S Ori 
25 suggests mass accretion. A surrounding accretion disc 
rotationally locked to the young central object could ex- 
plain the low angular velocity of S Ori 25. The angular 
velocity would be kept roughly constant at a low rate un- 
til the disc dissipates and the locking brakes. A similar 
picture is thought to occur in classical T Tauri stars. 

The periodograms of S Ori 27 and S Ori 28 display 
single peaks at 2.8 ± 0.4 and 3.3 ± 0.6 h, respectively. 
Folded WFC03 light curves clearly show sinusoidal modu- 
lations with amplitudes of about 30 mmag, as illustrated in 
Figure El Unfortunately, the WFC03 data cover roughly 
one cycle; hence, the determined periods are quite un- 
certain. The WFCOO light curves do not show periodical 
variations in time scales of ~ 3 h. We recall that S Ori 28 
is also a mid-term variable brown dwarf candidate. 

This kind of variability, possibly related to surface ro- 
tation, could be caused by heterogeneities spread over the 
photosphere (dust or spots) that evolve over timescales of 
several hours. The rapid evolution of atmospheric features 
that lead to the detection of "false" periods if only pieces 
of light curves are taken into account has been previously 
noticed by Bailer-Jones & Mundt (2001) and Gelino et 
al. (2002). Analogously to the Type Up classical T Tauri 
stars described by Herbst et al. (1994), hot spots linked to 
unsteady accretion or rotation could originate the modu- 
lations. Hot spots do not last for hundreds or thousands 
of rotations, as cool spots do. 

Other studies of photometric variability of a Orionis 
cluster substellar members have been published by Bailer- 
Jones & Mundt (2001) and ZO03. S Ori 31 and S Ori 45 
are the only objects in common with this work. Bailer- 
Jones & Mundt (2001) claimed a rotation period of 7.5 
± 0.6 h for S Ori 31. Although this timescale is suitably 
studied with the WFCOO data, the suggested amplitude 
of the photometric variation is quite small, ^lOmmag. 
These authors also found a peak of 30 ± 8 min in the pe- 
riodogram of S Ori 45, a ^25 Afj^p M8.5 brown dwarf with 
variable, strong Ha emission. A comparably short period 
(^46 min) has also been claimed by ZO03. However, as 
noted by those authors, given the expected radius of this 
young brown dwarf, such a short period cannot be asso- 
ciated with rotation because it would imply an extremely 
fast rotational velocity inconsistent with hydrostatic equi- 
librium. 

ZO03 also detected a second, more powerful, modula- 
tion in their optical and infrared light curves of S Ori 45 
with a period of 2.5-3.6 h that could be linked to rotation. 
The present WFCOO dataset has already been considered 
by ZO03. Our new analysis led to a very similar optical 
light curve. We retrieved the claimed period once the data 
points around phase 0.5 were excluded, as suggested in 
ZO03 after using the information from the infrared pho- 
tometry. According to our rather conservative criteria for 
the determination of photometric variability, S Ori 45 lies 
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Fig. 11. Amplitudes of variability versus I — J colour. 
Open triangles, squares and pentagons denote long-, mid- 
and short-term variables, respectively. Upper limits are 
shown for the rest of the objects in our sample. Strong 
Ha emission and near-infrared excess are also indicated. 



close to the borderline between variable and non- variable 
sources. We note that this is also the case of the periodic 
brown dwarfs S Ori 27 and S Ori 28. From our work and 
the data collected from the literature, we find that the fre- 
quency rate of periodic brown dwarfs in our final sample 
is around 18%. 

5.5. Final remarks 

Of our final list of twenty eight young brown dwarfs in 
the (J Orionis cluster, nine show some kind of optical pho- 
tometric variability (short-, mid-, long-term and periodic 
types) with a 5 % chance of a false positive detection, 
which represents 39 % of the sample. This frequency rate 
increases up to 46 % if S Ori 27, S Ori 28 (periodic short- 
term variables with slightly larger false positive detection 
probabilities), S Ori 31 and S Ori 45 (periodic brown 
dwarfs reported in the literature) are also considered. If 
we only take into account the a Orionis brown dwarf can- 
didates brighter than 7=19.5, for which the photomet- 
ric errors are reasonably low, the frequency of variability 
turns out to be 72 %. We summarize our results in TableEl 
where the abbreviations LT, MT and ST stand for long-, 
mid- and short-term variability, respectively. If a periodic 
signal is reported either in the present study or in the lit- 
erature, the letter "P" is used. S Ori 27, 28, 31 and 45 are 
considered as ST variables because of their short photo- 
metric periodicities, /-band amplitudes or scatters of the 
light curves and other interesting properties, like strong 



Ha emission and near-infrared excess, are also provided 
in Table ISl 

In Figure 1111 we plot maximum photometric ampli- 
tudes for the final sample of objects against the I ~ J 
colour, which is an indicator of temperature and mass. 
We use different symbols for the various types of vari- 
ability, and downward-point arrows for sources with no 
detected variability. As inferred from the figure, there is 
no clear trend between amplitudes and the colour of the 
brown dwarfs. S Ori J053825.4-024241 shows the largest 
amplitude, which, as discussed above, may be associated 
with the presence of a surrounding disc. It is also note- 
worthy from the figure that the largest /-band variabili- 
ties are found in the group of mid- and long-term variable 
brown dwarfs. We find /-band photometric variability in 
a Orionis cluster members with masses in the range 70- 
25Mj^,p. 

6. Conclusions and summary 

We have photometrically monitored in the /-band a sam- 
ple of thirty two young brown-dwarf candidates in the 
young star cluster a Orionis. Many of them are spec- 
troscopically confirmed as bona fide a Orionis substellar 
members. They are ~3 Myr old, and their masses span 
from the substellar limit down to the planetary-mass bor- 
derline. Twenty-eight of these substellar objects have re- 
liable variability indicators, and their light curves have 
been analysed in detail. This is one of the largest samples 
of young brown dwarfs ever monitored for photometric 
variability. 

With a 5 % chance of a false positive detection, we have 
found that nine of our targets show /-band variability on 
various timescales. If we include the periodic young brown 
dwarfs S Ori 27 and S Ori 28 found in our work, and S Ori 
31 and S Ori 45 reported in the literature, the incidence 
of variable u Orionis brown dwarfs is around 46 %. Such a 
frequency is higher than that (about 33 %) found among 
very low mass field stars and brown dwarfs. The measured 
amplitudes of variability are also more than three times 
larger in this study. This different behaviour can be as- 
cribed to the youth of the a Orionis brown dwarfs. We 
classify photometric variability on the basis of long-, mid- 
and short-term variations. 

— Long-term photometric variability (year-to-year varia- 
tions). Two substellar candidates (a minimum of 7 % 
of the final sample) show variations in the mean differ- 
ential magnitudes measured at two epochs 2.09 years 
apart, with amplitudes of 0.09-0.11 mag. They show 
neither inter- nor intra-night variations. 

— Mid-term variability (day-to-day variations). Six tar- 
gets (21 %) appear to be variable from night to 
night, with amplitudes of variation in the range 0.06- 
0.36 mag. Three of them have magnitude variations 
greater than 0.1 mag, which may make these the brown 
dwarfs with the largest amplitudes detected to date. S 
Ori 25 shows a ~40 h periodic modulation in the light 
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Name 




Variability 
type" 


Amplitude 

[mag]'^ 


Properties'^ 


S Ori J053847.2- 


-025756 


LT 


0.09±0.02 


- 


S Ori J053829.0- 


-024847 


- 


<0.03 


- 


S Ori J053954.3- 


-023719 


- 


<0.02 


- 


S Ori J053825.4- 


-024241 


MT + ST 


0.36±0.04 


nIR 


S Ori 25 




MT (P) 


0.15±0.02 


Ha 


S Ori J053826.1- 


-024041 


ST 


0.04 


- 


S Ori 27 




ST (P) 


0.03 


- 


S Ori J053922.2- 


-024552 


- 


<0.019 


- 


S Ori J054014.0- 


-023127 


MT 


0.12±0.05 


- 


S Ori 28 




MT + ST (P) 


0.07±0.03 


- 


S Ori 31 




ST (P'*) 


<0.02 


- 


S Ori 30 




- 


<0.015 


- 


S Ori 32 




MT 


0.06±0.02 


- 


S Ori J054004.5- 


-023642 


MT 


0.073±0.014 


- 


S Ori 36 




- 


<0.009 


- 


S Ori J053918.1- 


-025257 


ST 


0.05 


- 


S Ori 42 




LT 


0.11±0.03 


Ha 


S Ori 45 




ST (P'*) 


<0.03 


Ha 


S Ori J053929.4- 


-024636 


- 


<0.03 


- 


S Ori 71 




- 


<0.06 


Ha 


S Ori J053849.5- 


-024934 


- 


<0.05 


- 


S Ori 51 




- 


<0.04 


- 


S Ori 50 




- 


<0.05 


- 


S Ori 47 




- 


<0.04 


- 


S Ori J053944.5- 


-025959 


- 


<0.04 


- 


S Ori 53 




- 


<0.08 


- 


S Ori J054007.0- 


-023604 


- 


<0.07 


- 


S Ori J053956.8- 


-025315 


- 


<0.05 


- 



Variability type. LT: long-term; MT: mid-term; ST: short-term; P: periodic; -: no detected variability. 

Amplitudes of photometric variability: peak-to-peak amplitudes for LT and MT variables and ct(/) amplitudes for ST and 

non- variables. For the latter, the values must be understood as upper limits on variability. 

Properties or prominent spectral features. Ha: strong Ha emission; nIR: near-infrared excess. 

Periodic according to Bailer- Jones & Mundt (2001) and/or ZOOS. 



curve, consistent with the recent determination of the 
vsuYi parameter, which yields a rotational velocity of 
14 ± 4kms^^. Magnetic braking due to an accretion 
disc could account for the long period. 
Short-term variability (scale of hours). Three of our 
objects apparently display random non-periodic pho- 
tometric /-band variability on timescales of less than 
a few hours. The amplitudes of the light curves of two 
of these objects appear to be roughly steady during 
at least three consecutive nights with standard devia- 
tions at the 0.04-0.05 mag level. The third short-term 
variable brown dwarf, S Ori J053825.4-024241, is also 
a mid-term variable and shows the largest variability 
in our final sample. 

With slightly larger false positive detection probabili- 
ties, S Ori 27 and S Ori 28 display peaks in the WFC03 
periodograms at around 3 h, which could be related to 
fast rotation periods, such as those found S Ori 27 
and S Ori 28 additionally show peaks in the WFC03 



periodograms at around ^^3 h, which could be related 
to fast rotation periods, such as those found in other 
young brown dwarfs and evolved ultracool field dwarfs. 
These rapid rotations could be explained by the dis- 
sipation of a surrounding disc, or by the evolution of 
surface features on timescales of a few rotations. 

We note the correlation found between large ampli- 
tude of variability and the detection of a near-infrared 
excess or strong Ha emission, which could be an evidence 
of a substellar accretion disc. Four out of thirteen brown 
dwarfs with available spectra show strong Ha emission 
(i.e. pEW(Ha) > 30A). Two of them (S Ori 25 and S Ori 
42) are variable with amplitudes of ^^0. 12 mag and with 
timescales longer than one day. Regarding the other two, 
S Ori 45 has been independently confirmed as as photo- 
metrically variable by ZO03, and S Ori 71 is suspected to 
be a faint long-term variable. An infrared excess in the 
Ks band has been detected in the largest-amplitude vari- 
able brown dwarf of the sample, S Ori J053825.4-024241. 
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This is one of the first substellar objects sharing a corre- 
lation between near-infrared excess and large photometric 
variability, in analogy to low-mass T Tauri-like stars. 

A single variability scenario can hardly explain simul- 
taneously all the temporal scales involved. The presence of 
surrounding accretion discs, heterogeneous coverage over 
the photosphere by cool magnetic spots or fast-evolving 
dust clouds, mass transfer from or eclipses due to a com- 
panion have been suggested as possible causes of the vari- 
ability. Both long-term, multi-band coverage monitoring 
and complete spectroscopic studies are needed to conclude 
unambiguously which is the preferred variability scenario 
in young objects below the substellar limit. 
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